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Abstract 

The equilibrium magnetization for layered superconducting films that experience a nonzero 
component, H\\, of magnetic field applied parallel to the layers is computed at tempera- 
tures and at perpendicular field components in the vicinity of the decoupling transition. 
A fermion analogy is exploited for this purpose, whereby it is found that the parallel mag- 
netization shows an anomalous H^ 1 tail at high fields due to entropic fluctuations of the 
(parallel) lattice of Josephson vortices. A collective pinning effect is also identified for 
c-axis transport limited by a single planar defect oriented parallel to both c and to the 
applied magnetic field. 
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I. Introduction 



It is well known that the electronic conduction in the normal state of high-temperature 
superconductors (HTSC) is confined primarily to the copper-oxygen planes common to 
these materials. 1 This suggests the Lawrence-Doniach (LD) model of a stack of Josephson- 
coupled superconducting layers below T c . 2 ' 3 The following question can then be posed: 4 Is 
the superconducting phase strictly three dimensional (3D), or does two-dimensional (2D) 
superconductivity appear somewhere in the phase diagram? The experimental answer is 
complex. Conductivity measurements give evidence for a 2D superconducting transition, 5 
whereas thermodynamic measurements are more consistent with a 3D transition. 6 The 
resolution to this paradox is that, although the thermodynamic transition into the super- 
conducting phase is in fact 3D, 7-10 quasi-2D behavior arises for physics on length scales 
small in comparison to the Josephson penetration length, 11 Aj. Indeed, the c-axis Joseph- 
son plasma resonance seen in HTSC for magnetic fields oriented perpendicular to the layers 
that exceed the 2D-3D cross-over scale, 12 ~ ^o/A^, is an experimental example of such 
behavior, 13 in which case Josephson coupling persists in between layers while c-axis vortex 
lines have degenerated into decoupled pancake vortices. 3 

In this paper, we shall obtain the critical properties theoretically expected for thin 
films of extreme type-II layered superconductors in magnetic field aligned parallel to the 
layers. The quasi-2D regime cited above is presumed throughout. This guarantees that the 
thermodynamics "factorizes" into that associated with (perpendicular) pancake vortices 
and to that associated with (parallel) Josephson vortices. 9 ' 1 In addition, we shall assume 
that the thickness of the film (including leads) is less than the London penetration length to 
insure the absence of magnetic screening across the layers. The theoretical analysis exploits 
a new fermion analogy for the LD model. 15 We first compute the parallel equilibrium 
magnetization. It is found to exhibit a series of renormalized melting transitions connected 
with the parallel vortex lattice, 16 ' 17 as well as an anomalous H^ 1 tail at high parallel 
fields, 15 H\\ (see Figs. 1 and 2). We then compute c-axis transport limited by the presence 
of a single planar defect or pair of edges oriented parallel to both the c-axis and to the 
applied magnetic field. A collective pinning regime is identified in the high-field limit 
(see Fig. 1) for planar pins with an effective thickness £ p > a , where a denotes the 
lattice constant of the parallel vortex lattice. In particular, the current-voltage (I-V) 
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characteristic is found to be algebraic in this case, and to be inversely related to that 
obtained for a Luttinger liquid in the presence of a single backscattering impurity. 18 



II. Parallel Equilibrium Magnetization 



The problem to be solved then is the parallel thermodynamics of a finite number, N 
of Josephson coupled planes in the quasi-2D regime 11-14 H±_ > $o/Aj and/or T > T cr , 
where the 2D-3D cross-over temperature T cr marks the point at which the intra-plane 
phase coherence length, £, matches the Josephson penetration length, Xj. It is implicit 
then that long-range in-plane phase coherence is lost at a temperature T c below T* , at which 
point the planes entirely decouple; 7-9 i.e., £(T) = oo for T < T c , while Xj(T) = oo for 
T > T*. Since both £ and Xj are finite in the interval between T c and T*, the relationship 
£(T cr ) ~ Xj(T cr ) then implies the inequalities T c < T cr < T*. 11 We shall now introduce 
the corresponding LD free-energy functional, which reads 



in the absence of magnetic screening. Here, Q\ (r) denotes the phase of the superconducting 
order parameter in layer /, where r = (x, y) is the planar coordinate. The parallel magnetic 
induction B\\ = ($o/27rcZ)6|| found between layers / and l + l and aligned along the y axis is 



scale for the Josephson penetration length. The author has recently obtained as analogy 
between the above LD model and coupled chains of spinless fermions at zero temperature, 
where each chain corresponds to a layer. 15 Specifically, the Hamiltonian for the fermion 
model is divided into two parts, H = H\\ + H±, with 




(1) 



related to the vector potential above by A z = —b\\x. Here d represents the spacing between 
layers. Last, J\\ is a measure of the local in-plane phase rigidity, while A sets the bare 



H \\=J2 dx v F (^lid x ^ L -^ R id x ^ R ^+U\^l^ L ^ 



R 



(2a) 



1=1 



and 



H ± =U ± J2 / dx *L( x > l )*R( x > l + 1 )*L(x,l + l)* R (x,l) + 1i.c. , 




i=i 



(26) 
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and with field operators ^r(x, I) and ^l( x J) for right (R) and left (L) moving fermions. 
The coordinate along the Josephson vortices, y, is related to the imaginary time variable 
t of the fermion analogy by y = v' F r. Here, the Fermi velocity v' F = vf sech2</> is renor- 
malized by the intra-chain interaction t^, 19 with tanh2(/> = U\\/2ttvf- Also, U± > is 
a repulsive backscattering interaction energy 19 in between chains. The Gibbs free-energy 
of the LD model (1) with respect to the normal state is then found to be related to the 
ground-state energy E F of the fermion analogy by 

(G s - G n )/k B T = (L v /v' F )[E F (U ± ) - E F (0)). (3) 

The identifications 

b\\ =2n(N l+1 -N l )/L x , (4) 

T =e 2 ^ , (5) 

A " 2 =a- 2 {2\U 1 _\/nv' F ){T/T^), (6) 

complete the equivalence between the models, where Ni/L x gives the fermion density in 
the I th chain, T* = AttJ\\ is the decoupling temperature that marks the point at which 
interlayer phase coherence is lost, 7-9 and where uj g = i^a -1 is the ultraviolet cutoff in 
energy. This analogy is a direct generalization of the well-known equivalence that exists 
between the sine-Gordon model and the massive Thirring/Luther-Emery model in 1 + 1 
dimensions to a layered structure. 21 It reduces to a free theory in the double-layer case 
(N = 2) along the Luther-Emery (LE) line T = 2irJ\\. 19 

The general case, however, can be treated in the mean-field approximation defined 
by the charge- density wave (CDW) order parameter Xl( x ) = (^\i( x J)^L(xJ)), and the 
associated gap equation 22 

Ai = U\\xi + U ± (xi+i+Xi-i)- ( 7 ) 

[The order parameters at the boundaries are set to Xo( x ) = = xn+i(x)-] The mean-field 
Hamiltonian then has the form Hmf = J2iLi I dx^\{Hi — where the spinor field, 

^i(x) = {^>l{x, /), ^r{x, /)), for each layer is acted upon by the one-body operator 

H l =a 3 v F id x + a + A l (x)+a-A^(x). (8) 
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Here we define a± = \{(J\ ± «cr 2 ), where cr^ represent the Pauli matrices. Also, \i\ denotes 
the chemical potential for the fermions in layer /. At zero parallel magnetic induction, 
the mean-field ground state has spin-density wave (SDW) type order, xi( x ) = (~ l) Z Xo, 
for U± > 0. In the minimal case, A = 2, a constant gap Ao = u)o/smh[2nvF/(U± — U\\)] 
opens at each Fermi surface for interaction parameters satisfying U±_ > U\\. The latter 
phase boundary agrees with renormalization group results. 23 Also, by Eqs. (3) and (4), 
the line-tension of a single Josephson vortex is given in general by e\\ = kBT\A a \/v' F , 
with a pseudo-spin gap equal to A a = A (U±) — A (0). For weak coupling, U± — > 0, we 
then have A^ = U±g^-Ao\u ± =o at temperatures T < T* (or U\\ < 0), while A^ = A 
at temperatures T > T*o (or U\\ > 0). The former yields a pseudo-spin gap of A^ = 
0.9(U±/2na) along the LE line T = 2itJ\\ for the case N = 2, which is comparable to the 
exact value of A^ = U±/2wa. 19 The success that this mean-field approximation has in the 
double-layer case indicates that it is reliable for the case of large A, where fluctuations are 
smaller. A constant gap, A = cuo/sinh[27rfF/(2t/"_L — U\\)], opens at each Fermi surface 
in such case for interaction parameters satisfying U± > ^U\\. The identification (6) then 
implies the leading dependence T* = 47rJ||[l + |(a/A ) 2 ] of the decoupling temperature 
with the bare Josephson scale. Also, the pseudo-spin gap now has value A^ = 1.9(U±/2ira) 
along the LE line. The parallel lower critical field, = Ahe\\/<&o, vanishes exponentially 
at the decoupling temperature T#, however, since A a oc Aq(T). This implies an inflection 
point in its temperature profile below T*. 14,24 

The presence of a parallel magnetic field, however, will generally induce the phase 6[{x) 
of the CDW order parameter, xi = ( — l)'xo exp(z6' z / ), to wind. This can be seen explicitly 
from the Ginzburg-Landau (GL) equations that describe the present mean- field theory at 
"temperatures" T' = v' F /L y near T' c ~ A [see Appendix A, Eq. (A10)]. After extending 
Gorkov's original derivation 25 of the GL equations to the case of an isolated pair (/, / + 1) 
of consecutive layers, one recovers the original LD free energy functional (1), but with the 
bare scale Ao replaced by the renormalized Josephson penetration length Aj = vf /\A a \. 14 
This is a result of entropic wandering of the Josephson vortices in the parallel direction, 
which is particularly relevant in the regime of low parallel fields, b\\ <ti Aq -1 . Note that the 
later corresponds precisely to the long wave-length limit in which Gorkov's derivation of the 
GL equations is valid. 25 Also, the CDW phase is related to the true phase in the LD model 
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(1) by the gauge transformation 9[ = 9[ — 2kp,ix, where kp,i denotes the Fermi wave- vector 
of layer /. The formula Ni = ir~ 1 kF,iL x then indicates that the number of fermions, Ni, is 
equal to the winding number of the CDW phase in a given chain /. This fact coupled with 
the identification (4) also demonstrates that the gauge-invariant phase difference between 
these layers is just 9i + i — 9\ — A z = 9' 1+1 — 9[. A GL theory analysis 16 of the renormalized 
LD model 14 (1) then yields a sequence of first-order commensuration transitions at parallel 
fields h} between a parallel vortex lattice with flux penetration every l\ layers to one with 
flux penetration every /i + 1 layers as the field is lowered. Fig. 1 displays the predicted 
phase diagram in the critical regime using the estimates h\ = ^o/Xjd and h\ = 
based on the continuum limit. Contrary to claims in the literature, 8 we therefore find no 
evidence for a true layer decoupling transition as a function of parallel field in the critical 
regime. 14 ' 15 

Consider now a parallel vortex lattice in the limit of high parallel field, but where flux 
penetrates only every h > 2 layers. The LD model (1) indicates that spatial variations of 
the superconducting phase are generally absent in the former limit; 14 ' 16 i.e., 6i(x) = and 
9[{x) = 2kp,ix. If there exists no flux penetration in between layers / and / =F 1, the gap 
function (7) then has the form 

A z = (-l) l±1 U ±X oe i{e ' l±1 ~ e ' l) = (-l) l+1 A a e ±lb " x (9) 

along the special line U± = U\\ in parameter space, where b\\ denotes the flux in between 
layers / and / ± 1 and where = U±_xo- Here, we have made a gauge transformation in 
order to set ni and kp,i to zero (see ref. 22). The mean-field Hamiltonian (8) then has 
energy eigenvalues 

ef = v F k F ±[v 2 F (k-k F ) 2 + AI} 1 / 2 , (10) 

with an effective Fermi wave number kp = ±&y/2 and a pseudo-spin gap A a . The equilib- 
rium magnetization My = — gj^[(G s — G n ) /V] can then be computed using the equivalence 
(3). After following steps similar to those taken in the double-layer (/, / ± 1) case along the 
LE line, 15 one obtains the formula 



for the equilibrium parallel magnetization, where B\ = & /7t\jd is the parallel field cross- 
over scale, with Josephson penetration length Xj = vp/\A a \. Given the phase diagram 
(Fig. 1) arrived at by the previous low-field analysis, this result implies a series of jumps 
down in My at each melting field, Hf , as parallel field increases (see Fig. 2). Notice 
that the sign of the jump in this case is opposite to that observed in HTSC for vortex 
lattice melting as the perpendicular magnetic field is swept. 26 Last, (9) indicates that 
the gap function Aj (x) simply acquires a constant phase factor upon a uniform translation 
x — > x+ai of the coordinate in layer /. Hence, the parallel vortex lattice is infinitely smectic 
in the high-field limit. (This result is consistent with the shear instability obtained in the 
standard GL analysis of LD model (l). 3 ) The parallel magnetization in the high-field limit 
(I i = 1) is then given by that of an isolated double layer, 15 which coincides with (11) up 
to prefactor of 2 1 / 2 . This means that the parallel magnetization must show an anomalous 
H^ 1 tail at high fields. 

Finally, to demonstrate that the above formula (11) for the parallel magnetization 
in the high-field limit is generic, we shall now analyze the spectrum of the mean-field 
Hamiltonian (8) along another special line, U\\ =0, in which case degenerate perturbation 
theory in powers of U± can be employed. Let us suppose again that flux penetrates only 
every l± >2 layers, and that there exists no flux penetration in between layers / and / =F 1. 
After performing the same gauge transformation as before to set kp,i and \i\ to zero (see 
ref. 22), we obtain a gap function (7) of the form A; = (— l) z+1 A cr (l + e ±jb " x ), where 
by denotes the flux in between layers / and / ± 1, and where A^ = V j_xo- Now in the 
absence of Josephson coupling, U± = 0, the mean-field Hamiltonian (8) has eigenstates 
ty a = (0, e lkx ) and ^ b = (e^ fc_2fcF ^ x , 0), with the corresponding energy eigenvalues e a = 
vpk and Sb = vp^kp — k). As before, we have kp = ±b\\/2. The application of degenerate 
perturbation theory in powers of U± with respect to such states at momenta k ~ kp then 
yields the well-known formula 



4 = l( £ a + e b )±[^(e a -e b ) 2 + Al 



1/2 



for the perturbed energy eigenvalues. Substitution of the unperturbed energies above 
then yields the previous result (10). We thereby recover the formula (11) for the parallel 
magnetization in the case that h > 2. Last, when the parallel flux penetrates all layers 
(h = 1), the gap function (7) has the form A ; = A a (e^ lb w x + e ±lh w x ). This means 



that we must take into account the previous unperturbed states, as well as their time- 
reversed counterparts obtained after making the (global) replacement hp — > —k F . The 
end result is again the previous formula (11) for the parallel magnetization in the case 
h = l. 



III. Parallel Collective Pinning 

We shall now compute c-axis transport limited by a single planar defect or pair of 
edges in the high-field limit, H\\ > h\, where parallel flux penetrates in between every 
layer. The planar pin is assumed to be parallel to both the c-axis and to the applied 
magnetic field. Due to the extreme smecticity symptomatic of this regime, it is sufficient 
to study the double-layer case. The LD model (1) then reduces to a pinned sine-Gordon 
system 14 with free-energy functional 

E SG = J\\ J d 2 r[J(V0_ -£&||) 2 -Ao 2 cos#_] - e\ J dycos9.\ Xo 

+ elj d 2 r9(x-x )^, (12) 

where 9- is the gauge-invariant phase difference in between the consecutive layers. Here, 
(I/L y )($o/2ir) is the line tension due to the Lorentz force in the absence of flux 
flow. The Lorentz force and the pinning force are in equilibrium in such case, which implies 
that the c-axis current I flows only at the pin site xq. 27 Also, since Josephson coupling 
is weaker near the pin, we have that e p < 0. Following Coleman 21 and LE, 19 this model 
is equivalent to the massive Thirring model for ID fermion fields, ^> = (^l^r), m the 
presence of a single backscattering impurity. 18 Its Hamiltonian description is then 

H a = J dx^{a 3 v' F id x + a i A CT + 2 1 /V KF )* + 2g J dx^[^ R ^ L 

-£ P A CT (*I** + 4^)L, (13) 

where Vkf(x) = 2ivv' F (e^/kBT)9(x — x ) is the voltage drop equivalent to the Lorentz force 
located at xo, while £ p oc —€p gives the effective thickness of the pin plane. The interaction 
between fermions is related to the depinning temperature (or LE line) T^ p = 2nJ\\ by the 
relationship 

Zk = i + _^. (14) 

T 7TV' F K ' 



Notice then that the fermions interact repulsively for T < T dp , while they interact attrac- 
tively for T > T dp . Last, Eq. (4) indicates that the total number of fermions is equal to 
the total number of Josephson vortices lying in between the consecutive layers. 

Yet Eq. (14) indicates that the Thirring model is nearly free at temperatures near 
Tdp, with a Fermi surface at k F = \b\\. This suggests a canonical transformation of the 
form 

Ck = x k a k + y k b k (15) 
d k = -y k a k + x k b k (16) 

to remove the (mass) gap term in (13) at momenta k in the vicinity of the Fermi surface. 
Here, the original field operators for right and left moving spinless fermions are represented 
as ^r(x) = L x 1 ' 2 J2 k e lkx a k and ^l(x) = L% J2 k e lkx b k , respectively. The former is 
achieved by the choice of coherence factors (x k ,y k ) = (u k ,v k ) for k > and (x k ,y k ) = 
(v k , —u k ) for k < 0, with 

u k = 2 ^(1 + ^- 

, t = 2-^(l-^) V2 , (18) 

along with energy eigenvalue 



= 2-V(l + M)"', (17 ) 

)'pk 
E k 



E k = (v'*k 2 +Al) 1 / 2 . (19) 

Notice then that the energy eigenvalue for the new right-moving state state cjjO) is 
equal to e k = 9{k)E k — 9{—k)E k , while it is equal to — e k for the new left-moving state 
d\\0). Finally, the field operators obtained after such a canonical transformation are then 
^+(x) = Lx 1 ^ 2 J2 k e lkx c k and \l/_(x) = Lx 1 ^ 2 J2 k e lkx d k . After making the Luttinger 
liquid hypothesis, which assumes that only excitations near the Fermi surface are relevant, 
we obtain the following effective massless Thirring model Hamiltonian: 

H 'a = j dx^(a 3 v F id x +2 1/2 V KF )^' + 2g' J cfe#t#^# + #_ 

-Z' p A a (^_y + + y{y-)\ X07 (20) 

where v F = (v' F k F /E kF )v' F , g' = (v' F k F / E kF ) 2 g and £' p = (v' F k F /E kF )£ p are the 
renormalized Fermi velocity, interaction and pinning scale, respectively, and where = 



*+) is the canonically transformed spinor field. The above effective Luttinger model 
is valid in the limit of (i) a thick pinning plane, £ p > ao, with respect to the separation in 
between Josephson vortices, in the limit of (ii) high parallel fields B\\ > $o/A j<i, and at (Hi) 
temperatures in the vicinity of the depinning transition, T ~ Td p [see Appendix B, Eqs. 
(Bl) and (B14)]. Notice that the sheer existence of the gapless fermion analogy (20) indi- 
cates algebraic long-range order, (e^-^e -10 -^) oc r~ 2K °, for the vortex lattice, 20 where 
K' a = [2(T/T dp ) — l] 1 / 2 at high fields. Finally, transforming back the massless Thirring 
model (20) to the bosonic description, one recovers the gaussian limit of the sine-Gordon 
model (12) in zero parallel field, 

£ll = Jj J d\\{V9_) 2 -el J ^cos0_L + ^|4 J d 2 r9(x-x )^, (21) 

but with a renormalized local stiffness J|j such that 2nJ'^/kBT = 1 + g' /nv'p. Comparison 
of the latter with Eq. (14) thus yields 

^ = l + (22) 
k B T v F \T J y 1 

for one over the effective coupling constant of the gaussian theory (21), where 

-1/2 



n" 
V F 

V F 




(23) 



relates the Fermi velocity to the parallel field. 

Above, we have reduced the pinned sine-Gordon model (12) to the bosonic description 
of a Luttinger liquid with a single backscattering impurity. 18 In particular, the equivalence 
(3) yields the relationship Zsg oc Z^(U±) /Zll(0) between the respective partion func- 
tions. The linear density Xjl xn = (d y 9-/2ni) Xo for half-loop (or fluxon 7 ) excitations of 
Josephson vortices along the pin 28 is thus given by the difference [SlnZ^/S^ao)]]^ 1 - , where 
the pinning term is held fixed. Here, the field ao(y) is defined by d y ao = Vkf/v'f, where 
Vkf = 'Z'xv'pE^/kBT gives the voltage drop in the fermion analogy (13). Kane and Fischer 
have computed such functional derivatives using perturbation theory, 18 where they obtain 
the result Sin Z^/S^icio) oc (Vkf/Vo)^' with exponent [i' = Anjyk B T — l. Assuming that 
the voltage scale Vq in the fermion analogy is related to a (field dependent) current scale 



10 



Iq through the relationship V = 2nv' F eo/kBT, with So = c 1 (/o/-^y)(^ > o/27r), the above 
results lead to the relationship 





A«L« 77 " T ( 24 ) 



for the density of fluxons along the pin, with 

H = 2^-l (25) 

and Iq corresponding to the exponent and to the current scale, respectively, in the limit 
of infinite parallel field and/or Josephson penetration length. The assumption of flux 
creep dynamics 3 then implies that the true voltage drop between consecutive layers is 
V = (h/2e)(c/\fi xn ), where c is the average creep velocity. We therefore predict that 
the c-axis I-V characteristic limited by parallel collective pinning is algebraic at relatively 
high parallel fields. In particular, Eq. (24) yields a voltage drop V oc 1^ for low currents 
I — > 0. Also, both fj,' — > fj, and I' Q — > Iq as B\\ — > oo. Eq. (24) thus indicates that the 
system becomes more superconducting as parallel field increases. Finally, by Eq. (22) we 
have that fj, = 1 = fj,' at the depinning transition, T = Ta p . This means that ohmic flux 
flow sets in at temperatures above or equal to T dp by Eq. (24). Given the temperature 
dependence of the parallel field scale Bl ~ h\ (see Fig. 1), the c-axis current should then 
peak at some temperature T p < T dp for fixed field and voltage. 



IV. Discussion 



The previous results evoke the following image for the parallel thermodynamics of 
layered superconductors in the quasi-2D regime: Each Josephson vortex can be viewed as 
a string of width A confined to a given pair of consecutive layers. Parallel fluctuations 
of the string give rise to an effective Josephson penetration length, Aj(T) > Ao, as well 
as to entropic pressure. 14 ' 15 ' 24 The latter is responsible for both the anomalous -f^y -1 tail 
shown by the parallel magnetization (11) and for the parallel collective pinning effect 
discussed above. Although HTSC films that are equivalent to Josephson-coupled layers 
(1) already exist, their physical properties have been examined only at temperatures far 
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from the critical region. 17 ' 29 Comparable studies should be carried out within the critical 
regime of these materials to test the predictions made here. 

Finally, it must be stressed that all of the results obtained here are only valid for 
physics at length scales large in comparison to the ultraviolet cut-off a of the fermion 
analogy; e.g., for parallel fields B\\ < Qo/ad. Clearly the in-plane coherence length, £o, 
which is roughly equal to the size of a typical Cooper pair, provides a lower bound for a. In 
addition, the string image mentioned above suggests that the bare Josephson penetration 
length A supplies an upper bound for a. Where exactly within these limits a lies remains 
to be determined. 

It is a pleasure to thank S. Sorella, G. Gomez-Santos, H. Safar, M. Maley and U. 
Welp for very informative discussions. This work was supported in part by National 
Science Foundation grant DMR-9322427 and by the Spanish Ministry for Education and 
Culture. 
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Appendix A: Derivation of Ginzburg-Landau Equations for Coupled 
Charge-Density Waves in One Dimension 



Below, we shall recover the LD free-energy functional (1) from the mean- field ap- 
proximation (8) for the fermion analogy [Eqs. (2a) - (6)] by extending Gorkov's original 
derivation of the Ginzburg-Landau field equation 25 to the case of two coupled CDW states 
in one dimension. The limit of low parallel field, fry <C Aq 1 , is assumed. 

Consider two adjacent layers / = 1, 2 in isolation, each of length L y along the parallel 
field. Assuming periodic boundary conditions in this direction, we may then define the 
fictitious temperature T' = v' F /L y for the CDW state in the fermion analogy. Standard 
mean-field calculations then yield that the CDW is stable for fictitious temperatures below 
a critical temperature T' c = (e 7 /7r)A , where 7 denotes Euler's constant. Note that the 
thermodynamic limit L y — > 00 clearly corresponds to the fictitious temperature T' = 0. 
To begin the derivation, we first observe that the mean- field equation H^i = e^i plus 
the definitions 

G t (x, t; x', t') = - (T* R (x, /, t)V R (x', /, 0) (Al) 
Fi{x,t;x',t') =- (T^ R (x,l,t)^ L (x',l,t')) (A2) 

for the normal and the anomalous Greens functions, respectively, ultimately lead to the 
Gorkov equations 

(iuj n + iv F d x + ^i)Gi{x, x', iuj n ) + Ai(x)F*(x, x', iuj n ) = S(x - x'), (A3) 
(-iuj n - iv F d x + iii)F t *(x, x', iu n ) - A* l (x)Gi(x, x', iuj n ) = 0, (AA) 

along with the associated gap equations 

Al (x) = —U\\F*(x, t; x, t) — U ± F* (x,t;x,t), (Ah) 
A*(x) = -U\\FZ(x,t;x,t) - U±Ff(x,t;x,t). (A6) 

Here, t = iy/v' F is the fictitious time variable, while iu n denote the corresponding Matsub- 
ara frequencies. Applying Gorkov's method 25 to these set of equations yields the Ginzburg- 
Landau field equations 

-(d x - 2ik F +) 2 A\ + Vp 2 A 2 a (A\ + A*) = (A7) 
-(d x - 2ik F , 2 ) 2 A* + Vp 2 A 2 T (A{ +A*) = (A8) 
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to lowest order in A;, where 

4e 7 / 2ttv f U± \ . A . 

is the pseudo-spin gap. Above, ((z) denotes the zeta function. Finally, these equations 
can be integrated, the result of which is the Ginzburg-Landau free-energy functional 
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F0C I ^l^' 2 ^ 1 ) + l[^t~ 2kF ' 2 ) + ^7 2 [l - cos (<?i - ^)] I Ag, (A10) 

for gaps of the form Ai(x) = (—l) l A o e l0 ' l( - x \ Above, Xj = vf/\A<t\ is the Josephson 
penetration length. After making the gauge transformation 6[{x) = 0i(x) + 2kF,ix, we 
recover the form (1) of the LD free-energy. Last, it is worth mentioning that Aj and Aq 
are approximately equal along the LE line, T = 2ttJ\\, in the limit \U±\ \U\\\. This is 
demonstrated by observing that we have U\\ = —Gtivf/S and A = (Tvv' F /\U±\) 1 ^ 2 a in such 
case [see Eqs. (6)], and by substitution of the former into Eq. (A9). On the other hand, 
Eq. (A9) also indicates that Aj diverges exponentially as T approaches the decoupling 
temperature T* from below, since A = 2cuoexp[— 2hvf/(2U± — U\\)]. This agrees with 
results based on the Coulomb gas analogy 14 for the LD model (1), as well as with those 
based on a model for "frozen" layered superconductors in the Meissner phase. 24 
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Appendix B: Derivation of Effective Massless Thirring Model 

Below, we shall obtain the effective low-energy Hamiltonian (20) of the massive 
Thirring model analogy (13) for a pinned double-layer superconductor (12) in the presence 
of parallel field. The limit of a thick pinning plane 

Zp > ao, (Bl) 
with respect to the average separation in between Josephson vortices (fermions) is assumed 
to insure the validity of the Luttinger-Liquid hypothesis in the fermion analogy (13). Note 
that the above separation is related to the parallel magnetic induction and to the Fermi 
wavenumber by b\\ = 2-k / a$ = 2kp. The gapless Luttinger model will be achieved through 
the canonical transformation [see Eqs. (15) - (19)] of the original left and right moving 
fields 

9 R (x) = L" 1/2 £> fc c fc - y k d k )e ikx (B2) 
k 

V L (x) = L- 1 ' 2 J2(VkCk + x k d k )e ikx (S3) 
k 

into the new fields ^+(x) = L% 1 ^ 2 J2 k c ke lhx and = Lx ^ 2 J2 k d k e lkx . We now 

begin the derivation of Eq. (20) by applying this canonical transformation to each term 
in the fermion analogy (13). 

t f 

Kinetic Energy. The kinetic energy is given by H = J2 k (£ k c k c k + e^d' k d k ), where 
= ±9(k)E k ^f9(—k)E k are the energy eigenvalues of the quasi-particle excitations. Since 
only those excitations that are near the Fermi surface are relevant, we can approximate 
the quasi-particle energy spectrum by = Bp ± v'p(k =F kp), where v F = (y'pkp / E kp )v' F 
is the group velocity at the Fermi surface. This immediately yields the new expression 

H = J dxv%(*iid x V- - *\.id x *+) {BA) 

for the kinetic energy modulo a trivial shift of the chemical potential. 

Potential Energy. Using the form Vkf(^) = J2 q eiqx VKF{q) for potential energy drop 
at x = xo in terms of its Fourier transform Vkf(q), we can reexpress the corresponding 
term, Hkf = 2 1 / 2 J (IxVkf^^, i n the fermion analogy by 

H KF = 2 1/2 y~]v KF (q)[(xk+qXk + y k+q y k )cl +q c k - {x k+q y k - y k+q x k )c\ +q d k 



k,q 



{y k+q x k - x k+q y k )d\ +q c k + (y k+q y k + x k+q x k )dl +q d k ]. (B5) 
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After making the Luttinger-Liquid-type approximation (x k+q , Uk+q) —> (xkiUk) above for 
the coherence factors, which is valid in the limit (Bl) kp^ P 3> 1, we recover the original 
simple form 

H KF = 2 1/2 J dxV K F(y+y+ + (B6) 
for the potential energy in terms of the new fields. Here, we have used the identity 

4 + yl = 1 - 

Pinning/Backscattering. Let us first reexpress the pinning energy (13) as 

H pin = -A a / dx(*l* R + V r Vl). (B7) 



If we now take the long- wavelength limit (Bl) 3> 1, then the bounds on the above 

integral can be extended to ±oo. After Fourier transformation, we obtain the expression 

H pin -> -A a ^2[(x 2 k - y 2 k ){d\c k + c\d k ) + 2x k y k (c\c k - d\d k )] (B8) 
k 

as a result. If we then make the additional replacements 

x k -> u kp and y k -> (sgnk)v kp (B9) 

valid under the Luttinger-Liquid hypothesis, tracing back the previous steps leads to the 
effective low-energy Hamiltonian 

iy P in^-(<-<)e P A CT (*t* + +*t*-)Uo-2^^ F ^A CT (4*++^*-)Uo (sio) 

corresponding to the pinning term. Last, the second term above represents a trivial shift 
of the chemical potential in the long wave-length limit (Bl). The effective pinning term at 
the Fermi surface thus has the same form as the original, but with a renormalized pinning 
scale i' p = {v' F k F /E kF )i p . 

Interaction. Substituting the canonical transformation (B2) and (B3) into the forward 
scattering interaction term H\ = 2g J dx^ L ^ R ^R^L yields an equivalent expression of 
the form 

Hl = ~j~ (yk+qVkcl +q Ck + Vk+ q x k c\ +q d k + x k+q y k d\ +q c k + x k+q x k d\ +q d k )x 

X k,k',q 

x (x k i- q x k ,c\,_ q c k , - Xk'- q yk'cli- q dk' -yk'- q x k 'd\,_ q c kl + y k > - q y k >d\, _ q d k ,). 

(mi) 
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After making the replacements (B9) in the coherence factors above, which is valid for 
excitations near the Fermi surface, we recover the original form 

"FhA f ^T,t lT ,t 




for the interaction energy in terms of the new fields, in addition to a forward scattering 
contribution 

H 4 = 2g (u kF v kF f J2[P+(Q)P+(-q) + P-{q)P-{-q)] (513) 

q 

that appears in terms of the new particle-hole operators p+(q) = L x J2 k c\, +q Ck and 

— 1/2 t 

P-(q) = L x Ylkdk+qdk f° r right and left moving fermions. This interaction can be 
incorporated into the kinetic energy (B4) via Kronig's identity, 20 which yields the final 
result 



2 



n / ( v F k F 90 K \ . p , , , 

VF = V "\ a7 + 2^Bf-| {BU) 



for the effective Fermi velocity. Notice, however, that this correction is negligible in the 
limit of high parallel fields, <C E kp , and at temperatures near the depinning transition, 

1 5-0 | < TTVp. 
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Figure Caption 



Fig. 1. Shown is the phase diagram for thin-film layered superconductors in the quasi-2D 
regime as a function of the parallel field, H\\. The integer l\ designates a parallel 
vortex lattice in which flux penetrates every l\ layers (see ref. 16). Also, T c0 denotes 
the mean-field transition temperature of an isolated layer. 

Fig. 2. The equilibrium parallel magnetization [Eq. (11)] obtained from the mean-field theory 
approximation to the fermion analogy is displayed in the vicinity of the decoupling 
transition, T ~ T* and/or B± ~ $o/Aj. 
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